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ABREGE DESCRIPTIF " 

Precede de communications radiomobiles AMRT 
iteratif . 

Selon 1' invention, on repartit plusieurs sequences 
de reference dans chaque tranche temporelle allouee. A 
la reception, on a recours a une representation en 
treillis du canal multitrajets et a un precede 
iteratif . 

Application aux radiocommunicat ions avec les 
mobiles . 

25 Fig. 12 
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PROCEDE DE COMMUNICATIONS RADIOMOBILES AMRT ITERATIF 



DESCRIPTION 



Domaine technique 

La presente invention a pour objet un procede de 
communications radiomobiles AMRT ("acces multiple a 
repartition en temps", en anglais TDMA pour "Time 
Division Multiple Access") iteratif. 

Elle peut s'appliquer au systeme radiomobile 
europeen GSM [1] et son evolution GSM++ ou au systeme 
americain D-AMPS. 

Etat de la technique anterieure 

Les systemes AMRT permettent de repartir dans le 
temps plusieurs utilisateurs communiquant a travers un 
meme canal radioelectrique [1, 3] . On separe chaque 
utilisateur des autres en lui assignant periodiquement 
une tranche temporelle (ou intervalle, ou creneau 
temporel, ou "timeslot" en anglais) pour lui permettre 
d'emettre et une autre tranche temporelle pour lui 
permettre de recevoir. Pour se premunir des 
imperfections en terme de synchronisation temporelle 
entre emetteurs et de 1' interference entre tranches 
temporelles consecutives produite par la pluralite de 
trajets, des intervalles de garde sont prevus de part 
et d' autre de chaque tranche temporelle. 

La figure 1 annexee montre ainsi des tranches 
temporelles TAi, TA2, .y.. TAn allouees a N 
utilisateurs, avec des intervalles de garde IG. 
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Le signal radioelectrique emis durant une tranche 
temporelle est obtenu par transposition en frequence du 
signal en bande de base equivalent. Le signal en bande 
de base est le resultat du filtrage du train de donnees 
a emettre par un f litre dit d' emission. Habituellement , 
le train de donnees emis est compose de deux sequences 
de symboles de donnees separees en temps par une 
sequence de symboles connue du recepteur, dite sequence 
de reference. Les deux sequences de symboles de donnees 
peuvent eventuellement provenir du codage et de 
1' entrelacement des informations a emettre. 

Afin de permettre une estimation precise du canal 
au niveau du recepteur, la sequence de • reference 
possede souvent de bonnes proprietes de correlation 
temporelle. A titre d'exemple, les sequences de 
reference utilisees dans le cadre du systeme GSM, 
appelees sequences CAZAC (contraction de "Constant 
Amplitude - Zero Autocorrelation" en anglais) [1, 4, 
5] f sont des sequences a composantes prises dans 
1' alphabet bipolaire {-1, 1} et possedant une fonction 
d' autocorrelation circulaire nulle partout sauf a 
1' origine . 

Le canal radiomobile emprunte lors d'une 
communication entre un emetteur et un recepteur est 
generalement de type multi-tra j ets avec evanouissements 
rapides dits de RAYLEIGH. L' existence de plusieurs 
trajets vient de ce que I'onde radioelectrique se 
propage selon plusieurs chemins entre le lieu 
d' emission et le lieu de rece^ption, 

Le signal regu est alors la somme de plusieurs 
repliques plus ou moins retardees et plus ou moins 
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alterees en phase et en amplitude. Pour obtenir- une 
restitution fiable de 1' information transportee par ces 
repliques, le recepteur opere un filtrage adapte au 
filtre d' emission et au canal et il combine de maniere 
optimale les contributions energetiques de toutes les 
repliques du signal emis. Le signal en sortie du filtre 
adapte est echantillonne au rythme des symboles et il 
est blanchi au moyen d'un filtre discret dit 
blanchisseur . 

Les echantillons en sortie du filtre blanchisseur 
fournissent une version filtree et bruitee du train de 
donnees emis. Le filtre discret associe a cette version 
filtree, egalement appele canal discret, presente une 
reponse impulsionnelle finie qui varie d'un echantillon 
a 1' autre. II caracterise d'une maniere indirecte la 
realisation du canal mult i-tra j ets durant la tranche 
temporelle de reception correspondante . 

Les echantillons en sortie du filtre blanchisseur 
correspondant a la sequence de reference sont utilises 
pour estimer le canal discret [1, 4, 5]. Cette 
estimation du canal discret permet d'egaliser les 
echantillons restants permettant eventuellement de 
detecter les deux sequences de donnees emises. 
L'egaliseur habituellement utilise est connu sous le 
nom d'egaliseur de VITERBI avec decisions souples [2]. 
Comme son nom I'indique, cet egaliseur fournit, grace a 
une application de I'algorithme de VITERBI avec sorties 
ponderees ("Soft-Output Viterbi Algorithm" ou SOVA en 
anglais) [2], des decisior}s souples sur tous les 
symboles de donnees emis. Ces sorties souples sont 
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eventuellement desentrelacees et decodees pour detecter 
les informations envoyees. 

La figure 2 annexee illustre ces operations 
d ' emission/reception . L'emetteur E comprend une source 
de donnees 10, un codeur /entrelaceur 12 delivrant des 
symboles ak^ un modulateur 14. Le canal mult i-tra j ets 
est symbolise par le bloc 20. Le recepteur R comprend 
un demodulateur 30 (filtre adapte/f iltre blanchisseur ) 
delivrant des echantillons R^, un estimateur de canal 
discret 32, un egaliseur 34, un desentrelaceur /decodeur 
36 et finalement un destinataire 38. 

Le canal discret vu a la sortie du filtre 
blanchisseur peut varier de maniere significative d'une 
tranche temporelle a la suivante. Cette variation est 
due principalement au changement des conditions de 
propagation entre I'emetteur et le recepteur et a la 
stabilite en frequence au niveau du recepteur. 

Les conditions de propagation ont une influence 
directe sur le canal multi-t ra j ets observe. Elles 
changent en raison d'une modification de 
1' environnement ou d ' un deplacement de I'emetteur et/ou 
du recepteur. Elles entrainent une variation temporelle 
du canal discret, aussi bien entre tranches temporelles 
successives qu'a I'interieur d'une meme tranche. 

La variation du canal discret entre deux tranches 
temporelles successives allouees au meme utilisateur 
est d'autant plus grande que I'intervalle temporel 
entre ces tranches est grand. Cette variation est 
accentuee, meme au niveau d'une seule tranche 
temporelle, en raison d'une augmentation de la 
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frequence porteuse ou de la vitesse de I'emetteur ou/et 
du recepteur. 

Dans la pratique, la variation du canal discret 
entre deux tranches temporelles est suf f isamment grande 
pour empecher toute estimation adaptative de celui-ci. 
L' estimation de canal au niveau d' une tranche 
temporelle doit alors reposer sur les seuls 
echantillons de la sequence de reference 
correspondante . 

L' utilisation de cette sequence de reference 
fournit une estimation non biaisee du canal discret au 
milieu de la tranche temporelle. A I'instar du systeme 
GSM, on fait generalement I'hypothese d' un canal 
invariant au cours d' une tranche temporelle regue. Dans 
ce cas precis, 1' estimation obtenue au milieu de la 
trame peut servir sans degradation a 1 ' egalisat ion du 
reste de la trame. Mais, le deplacement a grande 
Vitesse de certains terminaux et la demande soutenue en 
services a hauts debits operant a des frequences 
radioelectriques de plus en plus grandes, rendent cette 
hypothese de moins en moins justifiee. Le canal discret 
peut en effet subir des variations signif icatives entre 
le debut et la fin d' une tranche temporelle donnee. 
L'ecart entre ce canal discret reel et son estimation 
se creuse de plus en plus lorsqu'on s'eloigne de la 
sequence de reference. Ceci peut entrainer une 
degradation importante et irreversible de la qualite de 
la reception et/ou des performances du systeme AMRT . 

Une solution [13] a ^deja ete proposee pour 
resoudre ce probleme de variation rapide du canal dans 
le cas des canaux radioelectriques lies aux 
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radiocommunicat ions par satellite. Cependant, cette 
solution ne traite que le cas d'un canal avec un seul 
trajet et suppose que les symboles, de reference ou de 
donnees, ont une energie transmise commune. Elle ne 
peut done pas etre appliquee aux canaux multi-tra jets 
rencontres dans les communications radiomobiles 
terrestres . 

Par ailleurs, a I'instar de 1' extension EDGE du 
systeme GSM, des debits importants peuvent etre 
obtenus, a bande f requentielle inchangee, par la mise 
en oeuvre de modulations utilisant un plus grand nombre 
d'etats. Le systeme devient alors tres sensible a la 
qualite de 1' estimation de canal, et il peut subir des 
degradations importantes meme a de tres faibles 
vitesses de deplacement, si la puissance et/ou la 
longueur de la sequence de reference ne sont pas 
augmentees . 

Le but de la presente invention justement est de 
remedier a ces inconvenients . 



Expose de 1 ' invention 

La presente invention a pour but d'accroitre les 
performances des systemes AMRT en reduisant le rapport 
entre la puissance regue et 1 ' interference d' acces 
multiple pour une qualite de reception donnee ou, en 
d'autres termes, d'ameliorer la qualite de reception a 
puissance de reception constante, Cette amelioration de 
la qualite permet, entre autres, d' augmenter la 
capacite et la couverture d^s systemes AMRT, tels que 
le systeme GSM et surtout son extension EDGE pour 
offrir des services a hauts debits. Cette amelioration 
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est obtenue notamment grace a une optimisation de 
1 ' estimation de canal au niveau du recepteur. Cette 
optimisation permet de contrecarrer toute degradation 
en performances engendree par une variation rapide du 
canal entre le debut et la fin d'une meme tranche 
temporelle. Elle permet egalement de reduire 
signif icativement las degradations engendrees par une 
augmentation du nombre d'etats de la modulation et ceci 
sans recours a une augmentation de la longueur ou de la 
puissance de la sequence de reference. 

Un autre but de 1' invention est de reduire le prix 
de revient des terminaux en les rendant, a qualite de 
reception inchangee, moins sensibles a 1 ' instabilite en 
frequence d' un oscillateur local bon marche permettant, 
entre autres, la transposition du signal regu en bande 
de base. 

L' invention permet encore de reduire, a qualite de 
reception egale, le nombre et/ou la puissance des 
symboles de reference transmis durant chaque tranche 
temporelle. Ce but est atteint par la prise en compte, 
dans 1' estimation du canal discret, des symboles de 
donnees qui sont en general plus nombreux que les 
symboles de reference. Ce but est egalement atteint par 
la prise en compte de maniere optimale des echantillons 
correspondant aux symboles de reference et/ou de 
donnees d'un nombre arbitraire de tranches temporelles 
allouees a 1 ' ut ilisateur . Cette prise en compte est 
effective meme en cas de variations importantes du 
canal entre deux tranches successives appartenant au 
meme utilisateur. Ce but est enfin atteint en prenant 
en compte de maniere optimale les echantillons 
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provenant du canal de signalisation commun. et 
d'eventuels utilisateurs sur la liaison descendants. 

Selon 1' invention, les symboles de reference sont 
repartis en groupes sur toute la duree de la tranche 
5 temporelle allouee, ce qui garantit une meilleure 
poursuite de la variation du canal et done une 
robustesse accrue . vis a vis des grandes vitesses de 
deplacement et de 1 ' instabilite en frequence de 
1' oscillateur local. Les sequences de reference peuvent 
10 etre choisies independamment les unes des autres, mais 
elles doivent toujours presenter de bonnes proprietes 
d' autocorrelation . Elles peuvent etre derivees de 
^i; sequences de type CAZAC de longueur reduite. 

]l\ Avec 1' invention, I'energie transmise par symbole 

I/I 15 peut varier d'un symbole a 1' autre sans que cela pose 

l^' des difficultes pour 1 ' est imat ion du canal, laquelle 

s'effectuera toujours de maniere optimale. 
i;i Selon 1' invention, le traitement d'une tranche 

temporelle se fait par bloc a chaque fois que les 
p 20 echantillons correspondant a cette tranche et 

eventuellement a ceux d'autres tranches (provenant du 
meme utilisateur, d'autres utilisateurs ou du canal de 
signalisation commun) sont disponibles pour effectuer 
1' estimation de canal. 
25 A I'instar d' un recepteur AMRT classique (tel que 

celui du systeme GSM) on commence toujours par une 
estimation grossiere du canal discret en utilisant 
uniquement des echantillons dependant des seuls 
symboles de reference de ^la tranche temporelle a 
30 traiter. Cette estimation grossiere permet de 
caracteriser , au niveau de chaque symbole de reference 



SP 16348. C/RS 



ou de donnee du bloc a traiter, symbole par symbole, 
1' evolution de la phase et de 1' amplitude de tous les 
coefficients (egalement appeles "trajets" par abus de 
langage) du canal discret retenus pour la realisation 
de 1' egalisation . 

L'une des caracteristiques de 1' invention est 
1 * utilisation d*un treillis pour representer le canal 
multitrajets . Cette caracteristique reduit la 
complexite du recepteur. En effet, si 1 ' on prend comme 
exemple des sequences de symboles de donnees comprenant 
n symboles^ sans codage, avec une modulation de phase a 
deux etats {MDP2) , chaque bit emis peut prendre 1 ' une 
ou 1' autre de deux valeurs. Comme ces valeurs sont 
equiprobables, cela signifie qu'il existe 2"" sequences 
de symboles de donnees possibles. La complexite d'un 
calcul traditionnel utilisant les probabilites 
partielles pour toutes les sequences possibles croit 
done exponentiellement (2'') avec le nombre de symboles 
par sequence. 

L * utilisation, conformement a 1* invention, d ' un 
treillis pour representer le canal reduit 
considerablement la complexite des calculs, car ceux-ci 
vont se subdiviser en sections, chaque section donnant 
lieu a un certain type de calcul, toujours le meme pour 
chaque section. La complexite augmente done comme le 
nombre de sections, lequel croit comme le nombre n. La 
complexite n'est done plus exponentielle mais lineaire. 

Les etiquettes des branches des sections du 
treillis changent au gre de , la variation temporelle du 
canal discret au niveau de chacune des tranches 
temporelles constitutives du bloc. Par ailleurs, un 
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algorithme different de 1 ' algor ithme de VITERBI .avec 
sorties souples (SOVA) [2] est utilise pour une 
caracterisation plus precise des probabilites des 
branches du treillis condit ionnellement aux 
5 echantillons regus . Dans le cas d'une mise en oeuvre 
optimale de I'estimateur de canal, un algorithme dit de 
BAHL ou BCJR (abreviation de "Bahl Cocke Jelinek 
Raviv") [6] peut etre utilise pour fournir ces 
probabilites conditionnelles . D'autres algorithmes [10] 
10 plus simples, provenant generalement d' une 

linearisation de cet algorithme, peuvent etre utilises, 

*fi avec cependant une legere degradation en performance. 

lis permettent de reduire considerablement la 

|il complexite de I'estimateur de canal, 

Ifl 15 Dans le cas d'un recepteur classique, les seules 

sorties ponderees fournies par 1' algorithme SOVA sont 
lS| directement exploitees apres desent relacement pour 

pi decoder les informations transmises. Dans le cas du 

i]f recepteur de 1' invention, les probabilites 

p 20 conditionnelles caracterisant indirectement les 
symboles de donnees sont exploitees, en plus des 
symboles de reference, pour fournir une estimation de 
meilleure qualite du canal discret. 

L' estimation amelioree du canal discret obtenue a 
25 la fin d'une iteration donnee peut etre utilisee pour 
augmenter la qualite et la fiabilite des probabilites 
conditionnelles des branches du treillis caracterisant 
le canal discret. Ces probabilites conditionnelles 
ameliorees peuvent, a leur tour, etre utilisees, 
30 . conjointement avec les symboles de reference, pour 
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apporter une amelioration supplementaire a 1 ' estimation 
du canal discret. 

En theorie, une estimation iterative ne conduit a 
une estimation optimale du canal discret qu'apres un 
nombre infini d ' iterations . En pratique cependant, 
quelques iterations suffisent sans grande degradation 
des performances par rapport a une estimation optimale. 

Une fois 1' estimation du canal discret accomplie, 
le treillis correspondant peut etre utilise soit par 
I'algorithme optimal de BAHL (ou tout autre 
simplification de celui-ci), soit par I'algorithme SOVA 
pour fournir des sorties ponderees ameliorees au 
decodeur . 

Le caractere optimal du recepteur selon 
1' invention est lie a la qualite de 1' estimation du 
canal discret vu durant une tranche temporelle donnee. 
Ce caractere optimal repose sur 1' utilisation d'un 
algorithme iteratif dit SAGE (abreviation de "Space- 
Alternating Generalized Expectation-Maximization 
Algorithm") [7] pour trouver la realisation de canal la 
plus vraisemblable conditionnellement aux echantillons 
regus du bloc a traiter. L' algorithme SAGE est une 
extension de 1' algorithme EM (contraction de 
"Expectation-Maximization") [8, 9] qui permet de se 
debarrasser du probleme de couplage entre coefficients 
du canal discret lors de leur estimation. 

L' estimation du canal discret repose egalement sur 
la decomposition de chaque trajet du canal discret a 
estimer selon un algorit^hme d' expansion dit de 
KARHUNEN-LOEVE [3]. Cette decomposition fournit, d' une 
part, une caracterisation souple des variations 
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temporelles de chacun des trajets discrets et, d' autre 
part, s'integre naturellement dans I'algorithme SAGE. 
Elle permet egalement de traiter le cas ou I'energie 
emise par symbole varie d' un symbole a 1' autre. 
5 L' estimation du canal discret repose enfin sur 

1' utilisation de I'algorithme de BAHL, qui fournit a 
chaque iteration de I'algorithme SAGE les probabilites 
des branches du treillis caracterisant le canal 
discret, condit ionnellement aux echantillons regus et a 
10 1' estimation de canal fournie a 1' iteration precedente. 

De fagon precise, 1' invention a done pour objet un 

i ?1 

•^=J precede de communications radiomobiles du type a acces 

Ul multiple a repartition en temps, dans lequel on alloue 

hi 

yl 15 a plusieurs utilisateurs des tranches temporelles pour 

II j 

1' emission et la reception de signaux radioelectriques , 
M et dans lequel : 

ill 

i^S - a 1* emission, on constitue, dans chaque tranche 

temporelle allouee, des symboles de donnees et 
20 des symboles de reference, on emet un signal 

radioelectrique contenant ces symboles de 
donnees et ces symboles de reference, ce signal 
empruntant un canal radiomobile mult i-t raj ets , 
- a la reception, on effectue un filtrage adapte 
25 du signal req;u, on echant illonne le signal 

filtre au rythme des symboles, on effectue une 
estimation du canal radiomobile, on traite les 
echantillons en tenant compte de cette 
estimation du canal et 1 * on restitue les 
30 symboles de donnees propres a la tranche 

allouee. 
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ce procede etant caracterise en ce que : 

- a 1* emission, pour chaque tranche temporelle 
allouee, on constitue plusieurs sequences de 
reference possedant des proprietes de correlation 

5 temporelle appropriees et 1 * on repartit ces 

sequences dans chaque tranche temporelle, 

- a la reception : 

i) pour estimer le canal radiomobile, on 
represente ce canal par un treillis 
10 comprenant des branches symbolisant les 

transitions entre deux etats successifs du 
canal et caracterisant les sequences 

IJi 

ill possibles a la sortie du canal, 

111 

m ii) on calcule les probabilites de ces branches 

j^J 15 en mettant en oeuvre un algorithme iteratif 

r: dans lequel on calcule, par un nombre fini 

^ (D+1) d' iterations et selon le critere de 

y probabilite maximum a posteriori, une 

I'U 

p: estimation d'une representation |G ^ ^j^^Q des 

20 L trajets de ce canal, 1* iteration de rang 

(d+1) (3d+l) permettant d'obtenir une 

reestimation {g^^'^'^^^ISTo ^ partir de 
1* estimation ^^^"^^f^'l obtenue par 1* iteration 
precedente (3d), et, a chaque iteration 
25 (3d+l), on calcule en L etapes les vecteurs 

jQ^(d+i)j^ ^=0, 1, L-1 a partir des 

vecteurs {G"^^^^»^}l"io et {c^^'^t'le calcules a 
I'etape precedente, 1' estimation {c^^^^j^J des 
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L trajets du canal etant finalement obtenue 
a partir de 1* estimation de sa 

presentation {g^^^^|^^q . 



re 



5 Breve description des dessins 

- la figure 1, deja decrite, montre la repartition 
des tranches temporelles allouees a plusieurs 
utilisateurs ; 

- la figure 2, deja decrite, est un schema 
10 illustrant 1' emission et la reception dans la 

technique AMRT ; 

^} - la figure 3 montre les variations du taux 

111 d'erreur binaire (TEB) en fonction de I'etalement 

A Doppler normalise (BdTs) ; 

15 - la figure 4 montre les variations du taux 

O d'erreur binaire (TEB) en fonction du rapport 

P| Eg/No pour des modulations iyiDP4 et MDP8, dans une 

premiere configuration et pour un canal discret a 

deux trajets ; 

20 - la figure 5 montre les variations du taux 

d'erreur binaire (TEB) en fonction du rapport 
Eg/No pour des modulations iyiDP4 et MDPS, dans une 
deuxieme configuration avec un canal discret a 
deux trajets ; 

25 - la figure 6- montre les variations du taux 

d'erreur binaire (TEB) en fonction du rapport 
Eg/No pour une modulation iyiDP4 utilisant la 
deuxieme configuration et un canal discret a 
quatre trajets ; 
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- la figure 7 montre les variations du -taux 
d'erreur binaire (TEB) en fonction du rapport 
Eg/N pour une modulation MDP4 utilisant un canal 
discret avec un, deux et quatre trajets ; 

- la figure 8 est une representation du canal 
discret vu a la sortie du filtre blanchisseur , 
sous forme de registre a decalage ; 

- la figure 9 illustre 1' estimation des 
coefficients du filtre du canal discret par 
correlation des echantillons regus avec une 
sequence CAZAC et ses decalees circulaires ; 

- la figure 10 montre la repartition des valeurs 
propres les plus importantes de la matrice de 
covariance pour trois etalements Doppler 
normalises ; 

- la figure 11 illustre une representation en 
treillis du canal discret normalise vu a la 
sortie du filtre blanchisseur pour un canal a 
deux trajets et une modulation a deux etats et 
une longueur de tranche temporelle egale a 
quatre ; 

- la figure 12 montre la chronologie des 
traitements effectues a la reception pour une 
tranche temporelle regue ; 

- la figure 13 montre la chronologie des 
traitements effectues dans une iteration 
elementaire de I'algorithme d* estimation de 
canal. / 
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Description de modes particuliers de mise en ceuvre • 

Les performances du procede selon 1' invention ont 
ete comparees a celles d'une variante du procede 
classique utilise dans le cadre du systeme GSM [l].Le 
recepteur classique du systeme GSM estime d' abord le 
canal discret au milieu de la tranche temporelle a 
traiter, par correlation des echantillons regus 
dependant uniquement de la sequence de reference avec 
la sequence CAZAC correspondante . II egalise ensuite 
les autres echantillons regus au moyen de 1 ' algorithme 
SOVA afin de decoder les symboles d' information emis. 
Afin d' evaluer equitablement 1 ' amelioration intrinseque 
de performances apportee par 1' invention par rapport au 
recepteur classique du GSM, I'algorithme sous-optimal 
SOVA a ete remplace par I'algorithme optimal de BAHL 
utilise dans 1' invention. 

La comparaison est obtenue a travers 1' evolution 
du taux d' erreur binaire (TEB) brut (sans prise en 
compte d' un eventuel codage correcteur ou detecteur 
d'erreurs), pour des symboles de reference aussi bien 
groupes au milieu de la tranche temporelle a traiter (a 
la maniere du systeme GSM) que repartis en petits 
groupes sur toute la duree de cette tranche 
(conf ormement a 1 ' invention) , afin de garantir une 
meilleure poursuite de la variation temporelle du 
canal. 

La figure 3 represente ainsi 1' evolution du TEB en 
fonction de I'etalement Doppler normalise B^T^, produit 
de I'etalement Doppler et de la duree d' un symbole 

, qui caracterise I'intensite de variation du canal 
due au deplacement de I'emetteur et/ou du recepteur. 
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L'etalement Doppler est a la fois proport ionnel "a la 
Vitesse du terminal mobile et a la frequence de la 
porteuse. L ' instabilite en frequence de 1' oscillateur 
local, qui n'est pas prise en compte ici, conduit a des 
degradations supplementaires en performances. A 
I'instar du systeme GSM, le comportement du TEB a ete 
evalue pour des tranches temporelles composees de N=148 
symboles chacune. La modulation retenue est la 
modulation MDP4 (modulation a deplacement de phase a 4 
etats) . Le canal discret considere est compose de deux 
trajets (L=2) de meme puissance moyenne. Le rapport 
Eg/M entre I'energie moyenne regue par symbole et le 
niveau de bruit est de 10 dB. Trois configurations ont 
ete considerees : 

■ Une sequence de reference centrale de 26 symboles 
bipolaires prise au hasard parmi les 8 sequences 
utilisees dans le systeme GSM. Ces sequences sont 
des extensions cycliques des 8 sequences CAZAC de 
longueur 16 [1, 11], 

■ Une sequence de reference centrale, {-1, +1, +1, 
+ 1, -1}, de 5 symboles bipolaires obtenue par 
extension cyclique (a un signe et un decalage 
pres) de 1' unique sequence CAZAC bipolaire de 
longueur 4 a savoir {-1, +1, +1, +1}. 

■ Cinq groupes de sequences de reference, repartis 
sur toute la tranche temporelle a traiter, 
identiques a la sequence de reference centrale de 
5 symboles bipolaires precedente (25 symboles de 
reference au total) , ' 
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recepteur classique, premiere configuration, 
recepteur classique, deuxieme configuration, 
recepteur de 1' invention, deuxieme 
configuration, 

: recepteur classique, troisieme 

configuration, 

: recepteur de 1' invention, troisieme 
configuration. 



Ces resultats ' montrent que 1 ' utilisation de 
sequences groupees au centre de chaque tranche 
temporelle (premiere et deuxieme configurations) ne 
permet d'aller au-dela d' un etalement Doppler normalise 
BdTs=1/2500 sans degradation - significative des 
performances. A 1' inverse, 1' utilisation de-groupes de 
petites sequences reparties (troisieme configuration) 
permet d'atteindre, sans degradation perceptible, des 
etalements Doppler normalises proches de BdTs=1/100. 
Cette amelioration significative de la robustesse du 
recepteur vis-a-vis de I'effet Doppler, et 
eventuellement de 1 ' instabilite en frequence de 
1 ' oscillateur local, est due a la repartition des 
sequences de reference qui permet de suivre de maniere 
plus precise toute evolution temporelle du canal. 

Le tableau 1 caracterise .de maniere chiffree 
1 ' amelioration en performances apportee par 1' invention 
vis-a-vis de I'algorithme cla^ssique ameliore du systeme 
GSM. Dans ce tableau R.C. signifie "recepteur 
classique" et R.I. recepteur selon 1* invention. Le 
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rapport Eg/NO entre I'energie moyenne regue par symbols 
et le niveau de bruit est de 10 dB. Dans le cas d'une 
sequence de reference centrale longue (premiere 
configuration), 1' amelioration apportee par 1' invention 
est negligeable. 



Etalement 


Premiere 


Deuxieme 


Troisifeme 


DOPPLER 


configuration 


configuration 


configuration 


normalise 


R.C. 


R.I. 


R.C. 


R.I. 


R.C. 


R.I. 
















1/2500 


0,00 
99 


0,00855 


0,01584 


0,0089 


0,0105 


0,0077 


1/100 


X 


X 


X 


X 


0,0143 


0,0083 



Tableau 1 



La sequence de reference, dans ce cas, est 
suffisamment longue pour garantir, sans I'aide des 
sequences de donnees inconnues, une bonne qualite 
d' estimation de canal. Cependant, dans le cas d'une 
seule sequence de reference courte (deuxieme 
configuration) , le recepteur selon 1' invention permet 
d'obtenir non seulement de meilleures . performances que 
le recepteur classique, mais aussi des performances 
comparables a ceux obtenus avec la sequence de 
reference longue de 26 symboles de reference. Par 
consequent, le recepteur selon 1' invention est rendu, 
d'une certaine maniere, insensible a une reduction meme 
importante de la taille de la sequence de reference. Ce 
resultat s'explique par la prise en compte, dans 
1' estimation de canal, des 'symboles des sequences de 
donnees qui viennent cooperer avec les symboles de 
reference. Dans le cas de cinq groupes de sequences de 
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reference courtes (troisieme configuration) , * le 
recepteur selon 1' invention permet de maintenir le taux 
d'erreur binaire presque inchange en faisant varier 
I'etalement Doppler normalise entre 1/2500 et 1/100. 

Les figures A, 5 et 6 representent le taux 
d'erreur binaire en fonction du rapport E^IN^ entre 
I'energie moyenne regue par symbole et le niveau du 
bruit pour un etalement Doppler normalise BdTs=1/2500 
peu contraignant . 

Pareillement au systeme GSM, le comportement du 
TEB a ete evalue pour des tranches temporelles 
composees de N=148 symboles chacune. Toutefois, les 
modulations retenues sont les modulations iyiDP4 et MDP8. 
Le canal discret considere est compose de un (L=l) ou 
de deux (L=2) ou de quatre (L=4) trajets de meme 
puissance moyenne. Deux configurations ont ete 
considerees : 

■ Une sequence de reference centrale de 26 symboles 
bipolaires prise parmi les 8 sequences utilisees 
dans le systeme GSM. Ces sequences sont des 
extensions cycliques des 8 sequences CAZAC de 
longueur 16 [1, 11]. 

■ Une sequence de reference centrale ({-1, +1, +1, 
+1, -1} pour deux trajets et {-1, +1, +1, +1, -1, 
+1, +1} pour quatre trajets de 4+L-l symboles 
bipolaires, obtenue par extension cyclique de 
1' unique (a un signe et un decalage pres) 
sequence CAZAC bipolaire de longueur 4 a savoir 
{-1, +1, +1, +1}. 
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Pour la figure 4, les modulations retenues sont- les 
modulations MDP4 et MDP8 utilisant la premiere 
configuration avec deux trajets (L=2) . Les courbes se 
rapportent aux cas suivants : 

- 51 : MDP8, recepteur classique, 

- 52 : MDP8, recepteur selon 1' invention, 

- 53 : MDP8, recepteur avec canal discret 

parf aitement connu, 

- 54 : MDP4, recepteur classique, 

- 55 : MDP4, recepteur selon 1' invention, 

- 56 : MDP4, recepteur avec canal discret 

parf aitement connu . 

Sur la figure 5, les modulations retenues sont les 
modulations MDP4 et MDP8 utilisant la deuxieme 
configuration et un canal discret avec deux trajets 
(L=2). Les courbes se rapportent aux cas suivants : 

- 61 : MDP8, recepteur classique, 

- 62 : MDP8, recepteur de 1* invention, 

- 63 : MDPS, recepteur avec canal discret 

par f aitement connu, 

- 64 : MDP4, recepteur classique, 

- 65 : MDP4, recepteur de 1* invention, 

- 66 : MDP4, recepteur avec canal discret 

parf aitement connu . 

Pour la figure 6, la modulation retenue est la 
modulation iyiDP4 avec la deuxieme configuration et un 
canal discret avec quatre trajets (L=4). Les courbes se 
rapportent aux cas suivants : 

1 
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- 71 : recepteur classique, 

- 72 : recepteur de 1* invention, 

- 73 : recepteur avec canal discret parf aitement 

connu . 

5 

Les courbes des figures 4 et 5 montrent que, a la 
difference d'un recepteur selon 1' invention, les 
performances du recepteur classique du systeme GSM sont 
tres sensibles a la longueur de la sequence de 
10 reference utilisee et se degradent substantiellement 
Q quand on passe de la premiere configuration a la 

1^ deuxieme. Les courbes des figures 5 et 6 montrent, en 

& outre, que les degradations en performances du 

y recepteur classique augmentent signif icat ivement avec 

15 le nombre de trajets alors que, a 1' inverse, ces 
degradations restent presque inchangees et faibles pour 
m un recepteur selon 1' invention. Le recepteur selon 

^ 1' invention presente done une excellente robustesse vis 

M a vis de toute reduction de la longueur de la sequence 

20 de reference et de toute augmentation du nombre de 
trajets du canal. 

Le tableau 2 caracterise de maniere chiffree, dans 
le cas de la modulation MDP4, 1 ' amelioration en 
25 performances apportee par le recepteur selon 
1' invention vis-a-vis de I'algorithme ameliore du 
systeme GSM. Le rapport BdTs est egal a 1/2500, le taux 
d'erreur binaire est egal a 10"^ et N=148. 
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Gain en dB 


Nombre de 
trajets ( £ ) 


Sequence de 
reference 
longue 


Sequence de 
r6f6rence 
courte 


Recepteur selon 1' invention avec 

sequence courte, recepteur 
classique avec sequence longue 


1 


0,03 


1,03 


0,03 


2 


0,48 


1,83 


0,48 


4 


1,42 


3,04 


1,43 



Tableau 2 



5 Ce tableau montre que, pour une longueur donnee 

(premiere ou deuxieme configuration) de la sequence de 
O reference, un recepteur selon 1' invention a toujours de 

't^! meilleures performances qu ' un recepteur classique 

III 

ameliore. Ces resultats montrent aussi que 

m 

y 10 1' amelioration apportee par 1' invention est d'autant 

^\ plus grande que la longueur de la sequence de reference 

est reduite ou que le nombre de trajets du canal est 
iii augmente. Ces resultats montrent enfin que, meme avec 

une sequence de reference courte (deuxieme 
O 15 configuration), un recepteur selon 1' invention a 

toujours de meilleures performances qu * un recepteur 

classique avec une sequence de reference longue 

(premiere configuration) . 

20 La figure 7 represente le taux d'erreur binaire en 

fonction du rapport E^IN^ pour un etalement Doppler 
normalise BdTs=1/100 contraignant et des tranches 
temporelles composees de N===100 symboles chacune. La 
modulation retenue est toujours la modulation MDP4 . Le 

25 canal discret considere est ^compose de un (L=l) ou de 
deux (L=2) ou de quatre (L=4) trajets de meme puissance 
moyenne. Une seule configuration avec trois sequences 
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de reference reparties est consideree. Chaque sequence 
est composee de 4+L-l symboles bipolaires et est 
obtenue par extension cyclique (a un signe et un 
decalage pres) de 1' unique sequence CAZAC bipolaire de 
longueur 4 a savoir {-1/ +1, +1, +1}. 

La valeur de I'etalement Doppler consideree est 
suf f isamment contraignante pour ne pas permettre a un 
recepteur classique de suivre les variations de chacun 
des trajets du canal discret qui peuvent subir des 
rotations de phase de ±7t entre le debut et la fin de 
chaque tranche temporelle. 

Les courbes de la figure 7 se rapportent aux cas 
suivants : 

- 81 : L=l, recepteur avec canal discret 

parf aitement connu, 

- 82 : L=l, recepteur de 1* invention, 

- 83 : L=2, recepteur avec canal discret 

parfaitement connu, 

- 84 : L=2, recepteur de 1* invention 

- 85 : L=4, recepteur avec canal discret 

parfaitement connu, 

- 86 : L=4, recepteur de 1' invention. 

Les courbes de ■ la figure 7 montrent que le 
recepteur selon 1' invention conduit a une degradation 
tres faible des performances vis-a-vis d' un recepteur 
ideal connaissant parfaitement le canal discret. Cette 
degradation est de I'ordre de 0,4 dB pour un seul 
trajet (L=l), 0,5 dB pour deux trajets (L=2) et 1 dB 
pour quatre trajets (L=4). 
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Les performances du precede de 1' invention ayant 
ete exposees, la phase de reception et notamment le 
traitement du signal regu, vont etre decrits plus en 
details . 

Un recepteur selon 1' invention determine, 
conditionnellement aux echantillons re?us du bloc a 
traiter, la realisation du canal discret la plus 
vraisemblable durant une tranche temporelle donnee. 
Cette determination se fait selon un critere dit de 
probabilite a posteriori maximale (MAP en abrege pour 
"Maximum A Posteriori probability" en anglais) . Cette 
estimation du canal necessite la connaissance des 
probabilites a priori des symboles des sequences de 
donnees transmises . 

Par souci de simplification, les developpement s 
mathematiques qui suivent et les figures 
correspondantes se rapportent toujours au cas ou chaque 
bloc a traiter pour 1' estimation de canal est compose 
uniquement des echantillons regus se rapportant a une 
seule tranche temporelle. Le cas ou les symboles (de 
reference ou non) d' autres tranches temporelles (du 
meme utilisateur ou d'autres utilisateurs , du canal 
commun de signalisation, . . . ) sont utilises se traite 
de la meme fagon. 

On note (Y I'operateur de transposition. Comme 
indique en figure 2, toute tranche temporelle d' un 
utilisateur donne est composee d'une suite de N 
symboles 3.^,3.^, que 1 * on peut ecrire sous forme de 
vecteur colonne 
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les symboles etant eventuellement modules en phase 
(MDP2, MDP4, MDPS, ...). L' energie transmise E^^ 
associee a chaque symbole a^^ est independante des 
valeurs prises par ce symbole mais peut varier d' un 
5 symbole de la tranche temporelle a 1' autre. Cette 
variation est generalement la consequence d' une 
allocation par le systeme AMRT d' energies differentes 
aux sequences de donnees et aux sequences de reference. 
Comme indique egalement dans la meme figure, les 
10 echantillons en sortie du filtre blanchisseur 
correspondant a cette tranche temporelle forment un 
I'l'i vecteur colonne 

W de N-vL-X composantes, ou L est le nombre de trajets 

111 15 du canal discret. 

La figure 8a rassemble les moyens compris entre 
^ 1' entree du modulateur de I'emetteur et la sortie du 

filtre blanchisseur du recepteur. Cette figure fournit 
done une representation, du canal radioel'ectrique reel 
20 et du traitement effectue au niveau du recepteur 
(demodulation, filtrage adapte et filtrage 

blanchissant) . Sur cette figure, on trouve le 
modulateur 91, un bloc 92 symbolisant le canal multi- 
trajets, le demodulateur 93 avec son filtre adapte 94 
25 et son filtre blanchisseur 95. 

La figure 8b montre une representation de ces 
moyens sous forme de registre a decalage avec (L-l) 
cellules Yi, Y2, Yl-i . Ce registre re?oit le symbole 

courant Sn et les cellules contiennent les symboles 
30 retardes d ' une periode symbole, a savoir an-ir 

an-(L-i)' Le circuit represents comprend encore (L) 



Q 
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multiplieurs Zq/ Zi, . . . , Zl-i recevant respect ivement 
(L) coefficients c^, cj,, ... c^"\ un addit ionneur S 
relie aux multiplieurs, un dernier additionneur 
recevant un signal Wn representant un bruit. Cet ultime 
5 additionneur delivre un signal rn- 

Chaque echantillon regu peut se mettre sous la 

forme 

/=o 

ou c[ est le /erne coefficient du f litre de canal 
10 discret au A: erne instant d' echantillonnage et Wj^ est un 
bruit additif blanc gaussien. Pour simplifier la 
III presentation et avoir des equations concises, les 

ill coefficients c[ du canal discret, verifiant k <l ou 

' k>N-\-l, interviennent dans 1' expression de R^^ alors 

p( 15 qu'ils sont nuls. Pour la meme raison, les symboles 

^\ emis, verifiant k<0 ou k>N-\, apparaissent dans 

cette expression de Rj^ meme s'ils sont nuls. 

Les coefficients du canal discret sont presque 
decorreles entre trajets discrets differents 
20 ( independance vis-a-vis de la variable /). Abstraction 
faite des effets de la modulation et de la 
demodulation, cette propriete est due a 1 ' independance 
des trajets multiples du canal radioelectrique reel. 
Toutefois, les coefficients d'un meme trajet discret / 
25 sont correles entre eux (dependance vis-a-vis de la 
variable k), Leur fonction d' autocorrelation discrete, 
definie par 

'pLn=E[cMy] 
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depend uniquement du spectre de puissance Doppler du 
trajet correspondant . 

Les spectres de puissance Doppler generalement 
rencontres en communications radiomobiles sont soit de 
type classique, (en environnements exterieurs) , soit de 
type plat, (en environnements interieurs) . La fonction 
d' autocorrelation du / eme trajet, de puissance moyenne 
(f>l , est donnee par 

<f>l=(l>oJ,{7tkB^T,), 

dans le cas d' un spectre de puissance Doppler classique 
et par 

(pl = (f>Q sm{KkB^T^) I TukB^T^ , 
dans le cas d'un spectre de puissance Doppler plat, ou 
Jq{') est la fonction de Bessel de premiere espece 

d'ordre 0 [3, 12]. La puissance moyenne varie d'un 

trajet a 1' autre et caracterise le profil de 
I'intensite des trajets ("multipath intensity profile" 
en anglais [3]) du canal discret. 

Le recepteur de 1' invention doit determiner selon 
le critere du maximum a posteriori (MAP) , la 
realisation la plus vraisemblable des vecteurs 
c' = (c' ,C/^, ,0,^2 9- • • ) ' / = 0,1,. . ., L — 1 , 

caracterisant les L trajets discrets retenus par le 
recepteur pour 1' egalisation . A cette fin, il a besoin 
d'une representation adequate du canal discret lui 
permettant de s'affranchir de la correlation temporelle 
des coefficients d'un meme trajet- II a egalement 
besoin d'une representat io/n qui lui permette de 
s'affranchir de la variation eventuelle de I'energie 
transmise d'un symbole de la tranche temporelle a 
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1' autre. Cette representation est obtenue tout d' abord 
par 1' introduction du vecteur de symboles transmis 
normalises 

avec des composantes Aj^ ^ I ^[E^ de module constant egal 
a un. Elle est egalement obtenue en reecrivant les 
echantillons regus durant une tranche temporelle sous 
la forme : 

/=o 

ou C[ est la A: erne composante du vecteur normalise 

associe et caracterisant directement le /erne trajet du 
canal discret. Elle est enfin acquise grace a une 
decomposition, dite de KARHUNEN-LOEVE [3], de chaque 

vecteur normalise C dans une base orthonormale |B j^^^ 
tenant compte a la fois .de la variation de I'energie 
des symboles emis et des caracteristiques du spectre de 
puissance Doppler du / eme trajet [3]. Cette 
decomposition permet en effet de decrire le vecteur 
normalise C' (avec composantes temporelles correlees) 
du /eme trajet discret par un vecteur colonne 
equivalent, G' = (Go,Gi',..,,Gjv.i)' , dont les composantes 
sont bien decorrelees. Le vecteur C' peut etre 
reconstruit a partir de sa representation ' adequate G' 
par la transformation 

*=0 
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Inversement , les composantes de cette 

representation adequate G' peuvent etre retrouvees a 
partir du meme vecteur C' par la transformation inverse 

Gl=(B''y'C\ /: = 0,1,...,A^-1, 

ou (•)* est I'operateur de conjugaison. 

Conf ormement aux proprietes de la decomposition de 
KARHUNEN-LOEVE [3], les vecteurs de la base 
orthonormale associee au / eme trajet ne sont autres que 
les vecteurs propres de la matrice de covariance H' du 
vecteur normalise associe a ce trajet. De pl.us, les 

valeurs propres associes ne sont autres que les 

variances des composantes Gq^G[,.,.,G^^_^ de la 
representation adequate G' du trajet correspondant . Les 
vecteurs {b'^}^_"J associes a chaque trajet, ainsi que les 
valeurs propres correspondantes , peuvent done etre 
obtenus a partir de 1' equations matricielle 

H'B'' =;L''B" , A: = 0,1,...,A^-1 . 

La (An,«)eme entree de la matrice de covariance H' 
du vecteur normalise associe au /eme trajet est donnee 
directement par 

ou 6^ denote la fonction d' autocorrelation discrete de 

T rn—n 

ce trajet definie plus haut. 

Le recepteur selon 1' invention doit normalement 
disposer de L bases orthonormales associees a chacun 
des L trajets du canal discr'et afin de reconstruire ou 
decomposer toute estimation intermediaire ou finale du 
canal discret. Dans la pratique, les formes des 
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spectres de puissance Doppler des L trajets du canal 
sont identiques. Les matrices de covariance {h'}^^_J 

associees sont done identiques, a un facteur 
multiplicatif pres, et les bases orthonormales 
correspondantes le . sont aussi. Une seule base 
orthonormale est alors generalement suffisante au 
niveau du recepteur pour effectuer convenablement 
1' estimation du canal. 

Pour generer la base orthonormale du /erne trajet, 
le recepteur doit disposer d'une estimation aussi 

precise que possible de la matrice de covariance H' de 
ce trajet. A cette fin, il doit disposer d' une 
connaissance aussi precise que possible des 
caracteristiques du spectre de puissance Doppler du 
canal radioelectrique . Cette connaissance est mise a 
jour regulierement et comporte obligatoirement une 
estimation de I'etalement Doppler Bd ou d'une borne 
superieure de celui-ci. Elle peut comporter 
eventuellement une estimation de la puissance moyenne 
de chacun des trajets du canal discret et de la forme 
du spectre de puissance Doppler. Conditionnellement a 
cette connaissance, le recepteur peut, soit calculer en 
temps reel les bases orthonormales adequates, soit 
rechercher dans une banque de bases orthonormales pre- 
calculees celles qui concordent le mieux avec les 
caracteristiques du canal reel. En general, le 
recepteur n'a ni connaissance precise de la forme du 
spectre de . puissance Doppler, ni connaissance parfaite 
et instantanee de I'etalement Doppler. Dans ce cas 
typique, il peut eventuellement faire I'hypothese d'un 
spectre de puissance Doppler plat avec, comme etalement 
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Doppler, une borne superieure de I'etalement Doppler 
reel. Par ailleurs, le recepteur peut ne pas disposer 
d'une estimation precise de la puissance moyenne de 
chacun des trajets du canal discret pris 
individuellement . Dans ce cas aussi, il peut aussi 
faire I'hypothese d'un canal discret avec L trajets de 
puissances moyennes egales. 

L' estimation des vecteurs c' caracterisant les L 
trajets du canal discret se fait iterativement selon le 

critere du MAP a travers les vecteurs normalises ou 
leur representation adequate G' . Cette estimation peut 
prendre en compte de maniere optimale aussi bien les 
caracteristiques des symboles de donnees et de 
reference emis ( stat istiques et energies) que les 
caracteristiques du canal multi-tra j ets . 

L' estimation au MAP • {g'}^^J d' une realisation de la 

representation adequate {G'}^^J du canal discret est la 
valeur 

(Gi:::=argmaj,4G'},';;iR) 

qui maximise la densite de probabilite condit ionnelle a 
posteriori Pi{^^}^~o\^^ ' Conformement a 1' invention, et 
en utilisant 1' extension dite SAGE [7] de I'algorithme 
EM [8, 9], il est possible d'obtenir iterativement une 
solution aussi proche qu' il le faut de la solution 
exacte . 

En general, la densite de probabilite 
conditionnelle a posteriori /'({G'}/^J|R) a maximiser 
possede plusieurs maxima globaux conduisant souvent a 
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une ambigulte de 1' estimation de canal selon le critere 
MAP. Cette ambiguite peut etre levee grace a 
1' utilisation des symboles de reference connus du 
recepteur. Cependant, ceci se revele souvent 
5 insuffisant, car cette densite de probabilite 
conditionnelle possede egalement des maxima locaux 
pouvant etre atteints par I'algorithme SAGE a la place 
de 1' unique maximum global. Pour eviter ce probleme, il 
faut utiliser les symboles de reference pour determiner 

10 convenablement les conditions initiales 

L'algorithme SAGE reestime successivement , par 
induction, chacun des vecteurs G' en gardant inchangees 

III les estimations precedentes des autres vecteurs C" , 

Ml 

111 avec m^l. Conformement aux proprietes de I'algorithme 

15 SAGE, ces reest imations successives garantissent une 

P croissance monotone de la densite de probabilite 

conditionnelle a posteriori p(|g'}^^q |R) . 

Etant donne le vecteur R d' echant illons regus 
durant une tranche temporelle a traiter, le recepteur 

20 commence par calculer des conditions initiales 

des vecteurs {Cj^^^Q a partir des echantillons ne 

dependant que des sequences de reference. A cette fin, 
les sequences de reference doivent etre telles qu'une 
separation rapide et simple des contributions des L 
25 trajets a chacun de ces echantillons soit garantie. 
Pour pouvoir effectuer une estimation initiale, 
permettant le lancement de 1'/ algorithme d' estimation de 
canal, il faut recourir a une ou plusieurs sequences de 
reference ayant chacune de bonnes proprietes 
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d' autocorrelation . A I'instar du systeme GSM, . ces 
sequences peuvent etre construites a partir de 
sequences CAZAC par extension ou permutation 
circulaire. Ces sequences sont composees de symboles de 
5 module constant et possedent une fonction 
d' autocorrelation circulaire nulle partout a 
1' exception de I'origine, Dans le cadre de la presente 
invention, il est preferable d' avoir des sequences 
CAZAC avec des symboles appartenant a la constellation 
10 de la modulation utilisee pour vehiculer les symboles 
de donnees. A titre d'exemple, pour une modulation 
iy[DP2, la seule sequence CAZAC (a un signe et a une 
permutation circulaire pres) la plus courte, a savoir 
U\ { -1, +1, +1, +1 } , ne compte que quatre symboles. Pour une 

tJl 15 modulation MDP4, utilisant la constellation {l,/ -1 ,-/} , la 

ill 

e; seule sequence CAZAC (toujours a un facteur 

K multiplicatif i, -1 ou et a un decalage circulaire 

|Jj pres) la plus courte, a savoir {l,/} / n'est composee que 

de 2 symboles. 

20 Une sequence de reference basee sur une sequence 

CAZAC de longueur M ne permet d'estimer que M trajets 
du canal discret au plus. Elle est obtenue par une 
extension cyclique de L-1 positions de cette sequence 
CAZAC et comporte au total iXI+L-l symboles. La longueur 

25 M de la sequence CAZAC doit done etre superieure ou 
egale au nombre de trajets L du canal discret a 
estimer. Dans le cadre de la presente invention, il est 
parfois souhaitable d'utiliser des sequences CAZAC 
courtes pour nombre de trajets L a estimer petit 

30 (inferieur ou egal a 4 pour une modulation MDP2 et a 2 
pour des modulations MDP4, MDP8, ...), car elles 
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permettent de reduire signif icat ivement le nombre total 
de symboles de reference par tranche temporelle, tout 
en garantissant un echantillonnage plus regulier et 
plus frequent du canal discret, surtout dans les 
environnement a grands etalements Doppler. 

La figure 9 illustre 1' estimation d'un canal 
discret compose de 2 trajets {L=2) dans le cas d'une 
modulation MDP2 . La premiere ligne, A, donne la 
sequence de reference, la seconde, Tl, correspond au 
premier trajet, la troisieme, T2, au second trajet et 
la quatrieme, S, donne la composition des echantillons 
regus. La sequence de reference elementaire utilisee 
{-1, +1, +1, +1, -1}, composee de 5 symboles (M+L-l = 5) 
n'est autre que 1' extension circulaire d'une seule 
position de la sequence CAZAC {-1, +1, +1, +1} composee 
de 4 symboles (iyi=4). Comme I'indique la figure 5, la 
sequence de reference elementaire occupe les positions 
entre n-4 et n de, la tranche temporelle a traiter. Les 
echantillons regus, dependant uniquement de cette 
sequence, occupent quant a eux les positions entre n-3 
et n. lis sont donnes explicitement par : 

R„., = C3 A-3 + Cl,A„., + IV„_, = C3 - Ci, + ^„.3 

R„., = C A„_, + Cl, A„_, + W„_, = C + d + IV„_, ■ 
R„ = C°A„ +Cl,A„ .+W„= -C° +C!,+W„ 

Pour une estimation peu biaisee des coefficients 
C° et Cl du canal discret normalise, pour k pris entre 
n-3 et n, il est necessaire / de faire I'hypothese d*un 
etalement Doppler suffisamment faible pour garantir une 
variation negligeable de canal entre les positions 
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et n. Dans ce cas, une estimation approximative de ces 
coefficients est obtenue en correlant les echantillons 
i?^ , k pris entre «-3 et n, avec la sequence CAZAC 
elle meme et un decalage circulaire d'une position vers 
5 la gauche de celle-ci. La correlation avec la sequence 
CAZAC {-1, +1, +1, +1} donne : 

comme premiere estimation des echantillons n-3,n-2,n-l 
et n du deuxieme trajet, ou est un terme qui tient 

10 compte du bruit et du biais. La correlation avec la 
sequence obtenue apres decalage circulaire, {+1, +1, 
+1, -1}/ donne quant a elle 

C°_3 = Q°_2 = = C'„ = (R„_, + R„_^ + R„_, -RJ/4 = 
(Cl, + C°_2 + C°_, + C°;/4 + 

comme premiere estimation des echantillons aux instants 

15 n-3,n-2,n-l et n du premier trajet, ou est un 

terme qui joue le meme role que . 

Dans le cas general, 1' estimation grossiere des L 
coefficients du canal, a un instant donne de la tranche 
temporelle, est obtenue par correlation des 

20 echantillons dependant uniquement de la sequence de 
reference avec L versions decalees de la sequence CAZAC 
correspondante . 

Chaque sequence de reference construite par 
extension circulaire d'une sequence CAZAC permet 

25 d'obtenir un echantillon bruite de chacun des trajets 
du canal discret. Le nombre 'd' echantillons par tranche 
temporelle necessaires a la reconstruction du canal 
depend de I'etalement Doppler normalise BqTs et du 
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nombre N de symboles composant cette tranche. Plus le 
produit BoTsN est grand, plus le canal varie d'une 
extremite a 1' autre de la tranche temporelle et plus le 
nombre de sequences de references elementaires doit 
etre grand. Une sequence de reference elementaire est 
un sous-ensemble de symboles groupes de la sequence de 
reference associee a la tranche temporelle obtenu a 
partir d'une extension circulaire d'une sequence CAZAC. 

Conformement aux proprietes de la decomposition de 
KARHUNEN-LOEVE du canal discret, plus une valeur propre 
A'* est grande pour un trajet / donne, plus la 
contribution G/b'* du vecteur propre associe a la 
representation adequate de ce trajet est importante en 
moyenne. La figure 10 montre la repartition des valeurs 
propres les plus importantes pour une tranche 
temporelle de N=100 symboles de meme energie normalisee 
a 1 et des etalements Doppler normalises BdTs egaux 
respectivement a 1/50 (courbe 113), 1/100 (courbe 112) 
et 1/500 (courbe 111). Cette figure montre clairement 
qu'une representation convenable de chacun des trajets 
du canal discret peut s'appuyer, sans perte perceptible 
de performance, sur les quelques vecteurs propres de 
plus grandes valeurs propres. Pour des valeurs faibles 
du produit BdTsN, le vecteur propre de plus grande 
valeur propre suffit a lui seul a une representation 
precise de chacun des trajets du canal discret. Ce 
vecteur propre varie alors tres lentement au cours du 
temps et permet done, grace a une sequence de reference 
groupee au milieu de chaque tranche temporelle, de 
reconstruire un canal discret quasi constant sur toute 
la tranche temporelle. Ceci est le cas du systeme GSM 
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ou les etalements Doppler normalises sont supposes etre 
tres faibles et le canal est suppose etre constant sur 
une tranche- temporelle donnee. 

Pour des etalements Doppler normalises un peu plus 
eleves, les deux vecteurs propres associes aux deux 
plus grandes valeurs propres sont necessaires a une 
reconstruction precise du canal discret. Ces deux 
vecteurs varient suffisamment pour qu'une sequence de 
reference composee de deux sequences elementaires , 
placees de part et d' autre du milieu de la tranche 
temporelle, se revele necessaire. 

Dans le cas general, plus grand est le produit 
BdTsN, plus grand est le nombre de vecteurs propres a 
utiliser dans la representation adequate du canal 
discret. De plus, plus grand est ce produit, plus 
importantes sont les variations de ces vecteurs propres 
et done plus important est le nombre de sequences de 
reference elementaires a mettre en oeuvre. De surcroit, 
a longueur de sequence de reference constante, plus le 
nombre de sequences de reference elementaires est 
grand, plus la longueur de chacune doit etre diminuee. 
Pour un nombre de trajets L donne, la longueur d'une 
sequence elementaire, M+L-1, peut au mieux etre reduite 
a 2L-1 en choisissant une sequence CAZAC de longueur 
L ; d'ou I'utilite d' avoir des sequences CAZAC tres 
courtes. Selon 1' invention, la sequence CAZAC {-1, +1, 
+1, +1} de longueur M=4 est recommandee pour toutes les 
modulations de phase quand le nombre de trajets L est 
inferieur ou egal a 4 et ^ que le produit BdTsN est 
faible. Selon 1' invention egalement, la sequence CAZAC 
{l,/} de longueur M=2 est conseillee pour un nombre de 
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trajets inferieur ou egal a 2 et des produits -BdTsN 
assez importants. 

Soit K le nombre de sequences de reference 
elementaires par tranche temporelle a traiter. Soient 
Cl^M^\^Cl_^^2r"^^k f / = 0,1,.,.,I-1 , les echantillons du canal 
discret procures par la A:eme sequence de reference 
elementaire et S,, = {k - M -hl.k - M -h2,,,.,k} I'ensemble des 
indices de ces echantillons. Les conditions initiales, 
|q/(0)||^-J^ de la representation adequate du canal selon 

1' invention peuvent etre obtenues a partir de ces 
echantillons en utilisant les equations : 

Le facteur de ponderation depend de la norme 

A^^ de la restriction des vecteurs propres aux seuls 
echantillons en correspondance directe avec les 
echantillons de canal procures par la sequence de 
reference. II depend egalement de la qualite de la 
contribution du vecteur propre B'* mesuree par le 
rapport ^!^INq. Plus la variance A'* du coefficient G[ 
dans la representation du canal discret est grande par 
rapport a la variance des echantillons de bruit, 

plus la ponderation correspondante est importante et 
s'approche de 1. 

En ' se basant toujours sur le vecteur R 
d' echantillons regus durant ^une tranche temporelle, le 
recepteur selon 1' invention conduit ensuite une 
estimation iterative du canal discret en partant de ces 
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conditions initiales {G'<°^},t;' . Conf ormement a 
1' extension SAGE de 1 ' algorithme EM, I'estimateur de 
canal optimal selon le critere du MAP calcule 
iterativement, trajet apres trajet, la reestimation 
j^i^-i ^ partir de 1' estimation {g'^''^},';' • A chaque 
iteration, I'estimateur de canal calcule de proche en 
proche les vecteur G'^''^') , / = 0,1,. z _ i , ^' partir des 
vecteurs {G"'<->}1-J„ et {g^^'^^'J, calcules auparavant. 
Selon 1' invention, la A:eme composante de la id + l)eme 
reestimation G'^''^'^ du vecteur G' associe au /erne trajet 
et obtenue en se servant de 1' expression : 

oil 

~ \ + No / A"' ' ~ ^^''-(i-»'^'.-(/.-2).---.'^J 



et C7''^ est la « erne composante de la reconstruction, 
C"'^'", selon la decomposition de KARHUNEN-LOEVE du weme 
trajet du canal discret a la d eme iteration. Cette 
reconstruction est donnee explicitement par 



% mk 



Comme precedemment , les ponderations wj. tiennent 
compte de la qualite de la contribution du vecteur 
a la representation adequate 'du canal discret. Plus la 
contribution d'un vecteur propre est importante par 
rapport a la variance du bruit, plus sa ponderation est 
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proche de 1 et plus il joue un role important dan-s la 
reconstruction de 1 ' est imat ion du canal. 

Afin d' estimer la representation adequate, G'^''''*^ 
du /erne trajet en accord avec 1' expression precedente, 
le recepteur commence par calculer une estimation 

m=0 /n=/ + l 

de C'v4 ; en retranchant les contributions des autres 
trajets m^l. Cette soustraction est faite dans 
I'hypothese ou le vecteur ^ (A„_^^_^^,A„-a-2)^-"^'^n) ^st 

emis, d' ou I'utilite de ponderer plus tard cette 
expression par la probabilite condit ionnelle a 

posteriori p(aJr,{G"^^"^^}|J^^4g"^^^};^^^ d' avoir reellement 

emis ce vecteur. Une fois cette estimation de C^A^^^ 
realisee, le recepteur la multiplie par le conjugue Al_i 
de A„_^ pour s'affranchir de la modulation. Cette 
operation est effectuee en supposant toujours que le 
vecteur A„ a ete emis, d' ou la necessite de ponderer 
1' estimation resultante : 

/-I L-\ \ 



fT = 0 /7l = /+l 

de Cl par la probabilite conditionnelle a posteriori 
d' avoir emis ce vecteur. L' estimation finale 

c:'"^" A„.„,- E c^"^ A„,„ k, /'(a„|r,{g'"-> } {G'"<''> ) ) 

de Cl est bien sur obtenue en prenant en compte la 
contribution de toutes les valeurs pouvant etre prises 
par A„ . 
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Une fois cette estimation , du « erne echantillon 
du /erne trajet obtenue pour n entre / est N + l-l, il 
ne reste plus au recepteur qu ' a la projeter sur chacun 
des vecteurs B* de la base orthonormale associee a la 
representation adequate de ce trajet. La projection sur 
le ^eme vecteur B* est ensuite ponderee par wl afin de 
tenir compte de la qualite de sa contribution a 
1' estimation du /erne trajet du canal a la id + l)eme 
iteration . 

En pratique, seuls quelques vecteurs de base 
interviennent de maniere importante dans la 
representation de chacun des trajets. Dans ce cas, le 
nombre de vecteurs utilises dans la transformation de 
KARHUNEN-LOEVE et la transformation inverse' peut etre 
reduit sensiblement sans perte perceptible de 
performances tout en garantissant une reduction 
importante de la complexite. 

L' estimation iterative du canal discret peut etre 
accomplie un nombre fini D de fois de telle maniere 

que 1' estimation globale {c'<°'}^^/J (obtenue a partir de 

{^'^°'}/=o 1^ transformation de KARHUNEN-LOEVE) 

garantisse une degradation imperceptible des 
performances du recepteur par rapport a la solution 
optimale qui serait obtenue apres un nombre infini 
d' iterations . Cette estimation iterative suppose la 
connaissance des probabilites conditionnelles a 

posteriori ^(a„|r,{G'"<''->}1":'^, {9'"<''>};;;;) pour toutes les 
valeurs possibles des vecteurs A„=(A„_^^_,^,A„.^,.2,,...,A„), 
n = 0,l,...,N + L-2 . Ces probabilites conditionnelles 
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peuvent etre calculees exactement par I'algorithme de 
BAHL, connu aussi sous le nom d' algor ithme BCJR [6]. 
Elles peuvent aussi etre calculees de maniere 
approximative et simple grace a des variantes 
linearises de I'algorithme de BAHL [10], 

A titre d' illustration, on presente ici les 
grandes lignes du calcul par I'algorithme de BAHL des 

probabilites condit ionnelles p(aJr,{g"^^"'^}|J^q,{G"*^*}^^^^^ . A 

cette fin, il faut remarquer, conformement a la figure 
8, que 1' echantillon i?„ obtenu a 1' instant n depend a 
la fois du symbole A„ emis a cet instant et des L-1 
symboles A^_^,A„_2,.,^,An.(L-\) 1^ precedant. Ainsi, 
1' ensemble modulateur-canal mult i-tra j ets-demodulateur 
peut etre modelise par une chalne de MARKOV a temps 
discret. Cette chaine de MARKOV est decrite par un 
treillis dont les etats representent a chaque instant n 
le contenu du registre a decalage de la figure 8 et 
dont les etiquettes des branches representent la 
version non bruitee des echantillons possibles a la 
sortie du f litre blanchisseur . A titre d'exemple, la 
figure 11 montre une representation en treillis d' un 
canal discret compose de 1 = 2 trajets pour une 
modulation MDP2 a deux etats {-1, +1} et une longueur 
de tranche temporelle N=4 . 

Un tel treillis represente-- la sortie attendue du 
canal en fonction du temps et selon les symboles emis. 
Dans le cas illustre ou le canal est suppose presenter 
deux trajets, la sortie du canal est conditionnee par 
le symbole emis a 1 * instant considere et a 1' instant 
anterieur. Par exemple, en partant de 1 * instant n=l, on 
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peut emettre soit -1, soit +1, sachant qu ' a 1 ' instant 
precedent 0, on n ' a rien emis (noeud initial a 1* extreme 
gauche). Si 1 ' on a emis -1, on s*attend en sortie du 

canal a recevoir (-1) Co+(0)Co , soit -Cq/ quantite qui 
5 figure sur la premiere branche montante vers le noeud 
-1- Si 1 * on a emis +1, on s' attend a recevoir 

( + 1) Cq+O Cq / soit Cq , quantite qui figure sur la 
premiere branche descendante vers le noeud +1. 

Pour n=l, si 1 ' on emet -1 et si 1 * on avait emis -1 
10 a 1* instant anterieur n=0, on s ' attend a recevoir 

(-1) C?+(-l)Cj , soit -C?-C} . Si I'on emet +1, on 

iff 

s'attend a recevoir ( + 1) C?"C}=c|'-ci . De meme, si 1 ' on 
Y\ avait emis +1 a 1' instant anterieur, on s'attend a 

J«J recevoir si 1 ' on emet -1, et cj^+cj si 1 * on emet 

15 +1. Ces quatre valeurs figurent sur les quatre branches 
p.i a n=l. La tranche s * acheve par une emission de 0, et la 

fij quantite attendue est -C4 ou C4 selon qu*on a emis -1 

pi ou +1 a 1* instant anterieur n=3. Le treillis represente 

done bien, a chaque instant, les quantites attendues en 
20 sortie du canal, selon les symboles envoyes. 

Conformement a une telle representation du canal 
discret, chaque branche du treillis determine une 
valeur unique du vecteur A„ = (^„_(^_,),y4^_(^_2),...,^„) / car, 
d'une part, les L-1 premieres composantes sont donnees 
25 par I'etat de depart S„ = (^„.(i.,jM,.(i_2)--- A,-t) de la 
branche et, d' autre part, les L-1 dernieres composantes 
sont donnees par 1/ etat de destination 

S„., =(^„-(L-2).^.-(i-3)---^J • Le probleme de calculer les 
probabilites condit ionnelles a posteriori 
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p(AjR,{G'"^^"'^}|^J^,{G'"^^^}^jj se ramene done au calcul des 

probabilites des branches du treillis 

conditionnellement au vecteur d' echant illons regus R 
et a 1' estimation intermediaire du canal discret 
normalise donnee par les vecteurs C'"^'''''^ pour les / 
premiers trajets w = 0,1,...,/- 1 et les vecteurs C"^^^ pour 
les L-£ derniers trajets w = /,/ + 1,. Z - 1 . A cette fin, 
les branches du treillis doivent etre etiquetees par 
les coefficients de cette estimation intermediaire et 
non par ceux du canal discret reel. 

L'algorithme de BAHL permet de calculer 
efficacement les probabilites conditionnelles a 
posteriori des branches du treillis en factorisant 
1' expression de la probabilite du vecteur regu 
conditionnellement a chacune des sequences A pouvant 
etre emises, afin de ne pas dupliquer toutes les 
operations communes a des sequences emises ayant le 
meme passe ou le meme avenir a un instant donne. A 
cette fin, cet algorithme fait appel a des fonctions 
a„(S) et y9„(S) qui donnent respect ivement un apergu du 
passe et de 1' avenir pour chaque instant n=0, 1, 
N+L-1 et chaque valeur S possible de I'etat S„ . 

La fonction a„(S) est la probabilite conjointe 
d'etre a I'etat S„=S a 1' instant n et d'avoir regu les 

echantillons - La factorisation des calculs pour 

des sequences ayant un passe commun permet de calculer 
cette fonction grace a la retpurrence : 

^„(S) = S^.-.(S'K(S',S) 



SP 16348. C/RS 



46 

7«(S',S) est la fonction qui donne la probabdlite, 
conditionnellement a un etat de depart du treillis 
S„., =S', d'aller a I'etat S„ = S et de recevoir en meme 
temps 1' echantillon . Cette expression montre que 

1' information sur le passe peut etre acquise a la fois 
a travers les probabilites associees a 1' instant 
precedant immediatement 1' instant actuel et a travers 
les probabilites de transitions des branches du 
treillis . 

La fonction fi„(S) a une signification equivalents a 
celle de a„(S) . C'est la probability d' avoir regu les 

echantillons |^4i=« ' conditionnellement a I'etat de 
depart S„ = S a 1' instant n. Cette fonction peut etre 
aussi calculee de proche en proche a partir de I'avenir 
par la recurrence 

A,(S) = S>9«.i(S)r„.,(s,s'). 

S' 

La encore, cette recurrence prouve que la 
connaissance de I'avenir pour un etat = S peut etre 
obtenue a partir des caracteristiques a la fois des 
etats S„^, =S' suivant immediatement cet etat et des 
branches du treillis les reliant ensemble. Pour une 
branche du treillis caracterisee par I'etat de depart 
S„ =S' = (^„_(^_i^,^^_(^_2j,...,^^^,) et I'etat d' arrivee 

S„^] = S = , la fonction ;Kjt(S',S) est donnee 

par : 



] f ^ /-I /-I 2^ 

X«.,(S',S) = i^(^J-— exp 

^^0 V ^^0 m=0 
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ou P{A„) est la probabilite a priori d'envoyer le 
symbole A„ . Ces probabilites a priori peuvent 
caracteriser un codage eventual des symboles modules 
reellement emis. lis peuvent egalement caracteriser les 
symboles de reference qui sont connus du recepteur et 
qui prennent une seule valeur possible. Pour ces 
symboles de reference, la probabilite a priori P{A„) 
prend la valeur 1 pour la valeur reellement emise et 0 
pour les autres valeurs possibles de la modulation. 

Afin d'effectuer les recurrences sur les fonctions 
a„(S) et A(S), il faut prendre comme valeurs 
initiales : 

ao(So =(0A-..,0)) = 1 et yff;v..-i(S^.,-i = (0A...,0)) = 1 

A partir du calcul par I'algorithme de BAHL des 
fonctions a„{S) et fi„iS) , il est facile d'obtenir les 
probabilites des branches du treillis, 

conditionnellement aux echantillons regus et a 
1' estimation intermediaire du canal discret, grace a la 
formule : 

«.(S')r..i(S',s)A.,(S) 



f(A. = A|R,{G-"-")i-:.,{G-"j:;;)= 



/?,(So=(OA.--,0)) 

ou A = (A-(/.-l)'^/i-(Z.-2)'---' A) / S' = A-(Z.-2)'---^«-l) 

II faut remarquer toutefois que la complexite de 
I'algorithme de BAHL reste toujours assez elevee et 
qu'en pratique il vaut mieux recourir a des algorithmes 
approximatif s qui permettent de diminuer notablement la 
complexite des calculs ^ans trop degrader les 
performances. Le recours a ces algorithmes 
approximatif s est d'autant plus conseille que le nombre 
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de trajets ou/et d' etats de la modulation . est 
important. 

Apres D iterations, le recepteur dispose d' une 
estimation du canal discret {C^^^^}^~^ proche de la valeur 
optimale selon le critere du maximum a posteriori qui 
est normalement obtenue apres un nombre infini 
d' iterations . A partir de cette estimation du canal et 
des • echantillons regus durant la tranche temporelle a 
traiter, le recepteur peut fournir soit des decisions 
fermes grace a I'algorithme de VITERBI, soit des 
decisions ponderees grace a I'algorithme SOVA ou a 
I'algorithme de BAHL. 

L'algorithme de VITERBI cherche la sequence emise la 
plus vraisemblable . II permet done de reduire la 
probabilite d' erreurs de decision sur ces sequences 
emises. L'algorithme SOVA est un algorithme de VITERBI 
avec sorties ponderees. II permet une caracterisation 
plus precise de la fiabilite des decisions effectuees 
par I'algorithme de VITERBI. L'algorithme de BAHL 
permet de determiner la densite de probabilite a 
posteriori de chaque symbole emis, condit ionnellement 
aux echantillons re?us et a 1' estimation de canal 
obtenue. II permet done de determiner la valeur la plus 
vraisemblable de chaque symbole emis et done de reduire 
le taux d'erreur. Les sorties ponderees sont obtenues a 
travers les probabilites des symboles 

conditionnellement aux echantillons regus. 

Dans tous les cas, les sorties ponderees ne sont 
utiles que si un decodage dorrecteur ou un codage de 
source avec entrees ponderees est utilise afin 
d'augmenter les performances du recepteur. 
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A titre d' illustration, dans le cadre. de 
I'algorithme de BAHL, le recepteur utilise 1' estimation 

canal discret pour calculer les fonctions 

a„(S) et /3„iS) . Ces fonctions permettent alors de 

determiner pour chaque symbole A„ sa distribution de 

probabilite par la formule : 

Toutes les operations precedentes d' estimation 
iterative du canal discret et de detection sont 
illustrees sur la figure 12. L'estimateur de canal 120 

utilise des moyens 121 definissant les bases [B ] , 
/ = 0,1,...,Z,-1 , et des moyens 122 calculant les 
coefficients de ponderation [w^ et \^jtj^=o * 

L'estimateur 120 met en oeuvre une premiere estimation 
grossiere du canal a partir des seules sequences de 
reference elementaires par correlation circulaire avec 
les sequences CAZAC correspondantes . Cette premiere 
estimation est symbolisee par le bloc 3o reference 
1201, qui delivre la representation adequate de 

1' estimation initiale {g'^^^}^^q • Ensuite, 1' estimateur 
120 met en oeuvre une reestimation d' ordre 
representee par le bloc 3^ reference 1202, qui delivre 
la representation adequate de 1' estimation de canal 

{G'^^*}^^J"J , suivie d' une reestimation d' ordre d + \, 
representee par le bloc 3^,i reference 1203, qui delivre 
la representation de 1' estimation de canal {G'^''''**},^q . 
Enfin, I'estimateur 120 met en oeuvre une derniere 
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reestimation d' ordre D, representee par le bloc 
reference 1204^ qui delivre la representation adequate 

de 1' estimation de canal {G'^^^}/to' • ^ partir de cette 
representation adequate du canal estime, 1 ' est imateur 
fournit par les moyens 1205 1' estimation definitive 

{C'^^^]^~Q du canal discret en mettant en oeuvre la 

transformation de KARHUNEN-LOEVE . 

La fourniture a un decodeur correcteur d'erreur ou 
a un codeur de source de decisions fermes ou souples 
sur les symboles emis, s^effectue par des moyens 123 
qui mettent eventuellement en oeuvre I'algorithme de 
VITERBI [3], I'algorithme SOVA [2], I'algorithme de 
BAHL [6] ou I'algorithme de BAHL linearise [10]. 

Selon 1' invention, la reestimation iterative du 
canal est effectuee D + 1 fois et permet de calculer 

' {G'^'^lt;' . ■•■r {G'^^^lt'o' , successivement . La 
figure 13 illustre les L etapes composant chaque 
iteration et qui sont necessaires a 1' estimation 
successive de chacun des L trajets. Grace aux moyens 
121 et 122 de la figure 12, (fournissant les bases 

{^'"^}r=o* ponderations {w^j^jj ) , la (c/ + 1) eme 

iteration 3ci+i I'estimateur de canal met en oeuvre 

une premiere etape, basee sur I'algorithme de BAHL et 
representee par le bloc BAHLq reference 130, pour 
delivrer les probabilites conditionnelles a posteriori 

p|a„|r,{G"'^''^)^^^'|jJ a partir de la d erne reestimation du 

canal. A partir de ces probabilites conditionnelles, la 
(irf + l)eme iteration met ensuite en oeuvre une deuxieme 
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etape, basee sur 1' algorithme SAGE et representee- par 
le bloc SAGEq reference 131, pour permettre la (d + \)eme 

reestimation G°^''^^* du premier trajet. 

A partir de la (<^ + l)eme reestimation des /-I 
5 premiers trajets, la (d-hl)eme iteration met en oeuvre 
une troisieme etape, basee de nouveau sur I'algorithme 
de BAHL et representee par le bloc BAHL/ reference 132, 
pour delivrer les probabilites condit ionnelles a 

posteriori p[^n\^lG^^'^'^y^^^^^^^ ^ partir de la d erne 

O 10 reestimation des L-l derniers trajets du canal et de 

la (c/ + l)eme reestimation des / premiers trajets de ce 
canal. Par la suite, la (d + l)eme iteration met en 

J^f oeuvre une quatrieme etape, basee sur I'algorithme SAGE 

III 

Rl et representee par le bloc SAGE/ reference 133, pour 

CJ 15 permettre la (<i-t-l)eme reestimation g'^"^""*^ du /erne 

d trajet. 

Enfin, a partir de la (^f + l)eme reestimation des 
^-•^ L-\ premiers trajets, la (d'hl)eme iteration met en 

oeuvre une derniere etape, basee sur I'algorithme de 
20 BAHL et representee par le bloc BAHL^_, reference 134, 
pour delivrer les probabilites condit ionnelles a 

posteriori p(aJr,{G'"^''*^^};^;J^,G^'-'^'^^^) a partir de la ^ eme 

reestimation du dernier trajet du canal et de la 
((3f-fl)eme reestimation des L-\ premiers trajets de ce 
25 canal. Par la suite, la {d-hl)eme iteration met en 
oeuvre une derniere etape, basee sur I'algorithme SAGE 
et representee par le bloc SAGE^_, reference 135, pour 
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permettre la (i/ + 1) erne reest imat ion g*^"*^"^*^'''^ du -Leme 
tra j et . 

La premiere iteration, symbolisee par le bloc 3o 
reference 1201 de la figure 12, ne peut pas mettre en 
oeuvre les algorithmes de BAHL et SAGE pour calculer la 
premiere estimation du canal discret. Selon 
1' invention, on utilise alors les symboles de reference 
dont les probabilites conditionnelles a priori et a 
posteriori sont bien connues du recepteur pour 
effectuer cette premiere estimation. Cette premiere 
estimation s'effectue, en correlant les echantillons 
regus, dependant uniquement de ces symboles de 
reference, avec les sequences CAZAC correspondantes . 
Elle met en oeuvre les moyens 122 de la figure 12 pour 



calculer les ponderations 
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REVEND I CAT I ONS 



1, Procede de communications radiomobiles du type 
a acces multiple a repartition en temps, dans lequel on 
alloue a plusieurs utilisateurs des tranches 
temporelles pour 1' emission et la reception de signaux 
radioelectriques , et dans lequel : 

- a 1* emission, on constitue, dans chaque tranche 
temporelle allouee, des symboles de donnees et 
des symboles de reference, on emet un signal 
radioelectrique contenant ces symboles de 
donnees et ces symboles de reference, ce signal 
empruntant un canal radiomobile multi-tra j ets , 

- a la reception, on effectue un filtrage adapte 
du signal regu, on echantillonne le signal 
filtre au rythme des symboles, on effectue une 
estimation du canal radiomobile, on traite les 
echantillons en tenant compte de cette 
estimation du canal et 1 * on restitue les 
symboles de donnees propres a la tranche 
allouee, 

ce procede etant caracterise en ce que : 

- a 1' emission, pour chaque tranche temporelle 
allouee, on constitue plusieurs sequences de 
reference possedant des proprietes de correlation 
temporelle appropriees et 1 ' on repartit ces 
sequences dans chaque tranche temporelle, 

- a la reception : 

i) pour estimer le^ canal radiomobile, on 
represente ce canal par un treillis 
comprenant des branches symbolisant les 
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transitions entre deux etats successifs du 
canal et caracterisant les sequences 
possibles k la sortie du canal, 
ii) on calcule les probabilites de ces branches 
en mettant en oeuvre un algorithme -iteratif 
dans lequel on calcule, par un nombre fini 
(D+1) d' iterations et selon le critere de 
probabilite luaxdLmum a posteriori, une 

estimation d'une representation 
L trajets de ce canal, 1' iteration de rang 
(d+1) (3d+l) permettant d'obtenir une 

reestimation {G'^^'^^^j^lJ ^ partir de 
1' estimation ^^^"^^f^^l obtenue par 1 ' iteration 
precedente {3d), et, a chaque iteration 
(3d+l), on calcule en L etapes les vecteurs 
QfCd+i)^ ^=0, 1, L-1 a partir des 

vecteurs ^'"^'^^^I'io {c'^^^^t'j, calcules a 

i*6tape precedente, 1* estimation {c*^^^}^J des 
L trajets du canal etant finalement obtenue 
a partir de 1' estimation de sa 

representation {c^^^^j^^o . 

2- Precede selon la revendication 1, dans lequel, 
a chaque etape de chaque iteration, on met en oBuvre 
successivement : 

i)un algorithme de BAHL ou une de ses variantes 
pour deiivrer des probabilites 

conditionnelles a partir de la a 
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reestimation des L-^ dernier s trajets et de 
ia (d+1)^^"*^ reestimation des i premiers 
trajetSr 

ii) un algorithrae SAGE pour delivrer une (d+l)"-^ 
reestimation du i^^"^ trajet. 

3- Precede selon la revendication 1^ dans lequel;- 

ayant obtenu une estimation du canal discret |c^^^^}^^o 
et, compte tenu des echantillons re<?us durant la 
tranche temporelle traitee, on prend une decision aur 
la sequence des symboles emis en utilisant un 
algorithme optimal de BAHL, ou une de ses variantes, 
pour obtenir la probabilite a posteriori de chaque 
symbole emis . 

4. ProcedS selon la revendication 1, dans leqiael 
on r^partit K sequences de reference dans chaque 
tranche allouee avec K compris entre 2 (inclus) et 12 
( inclus) , 

5. Precede selon la revendication 1, dans leqael 
chaque sequence comprend un nombre de symboles de 
reference inferieur 4 8. 

6- Precede selon 1 ' une quelconque des 
revendications 1 a 5^ dans lequel les sequences de 
reference possedent une fonction d ' autocorrelation 
circulaire nulle partout sauf a I'origine. 

7. Precede selon la revendication 6^ dans leqael 
on utilise une sequence de reference {-If +lf +lf +1} 
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ou une sequence derives de celle-ci par permutation 
circulairer ou changement de signe pour un nombre de 
trajets (L) inferieur ou egal a 4. 

8- Procedfe selon la revendication 6, dans leqviel 
on utilise une sequence de reference {1/i} ou une 
sequence derivee de celle-ci par permutation 
circulaire, conjuguaison ou changement de signe, pour 
un nombre de trajets (L) infferieur ou fegal a 2. 
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